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Abstract . We show that Alfven waves or Alfvenic surface waves can carry 
enough energy into the corona to provide the coronal energy requirements. 
Coronal loop resonances are an appealing means by which large energy fluxes 
can enter active region loops. The wave dissipation mechanism still needs to 
be elucidated, but a Kolmogoroff turbulent cascade is fully consistent with the 
heating requirements in coronal holes and active region loops. 

Introduction 


The solar chromosphere and corona are heated mechanically. The energy 
requirements of the chromosphere and corona are roughly comparable but one 
usually speaks of 'the coronal heating problem* , presumably because of the 
spectacularly high temperatures there; it is probably a mistake to conceptu- 
ally separate the chromospheric and coronal heating problems, but space 
constraints require us to do so here. Spicules present the dual problems of 
heating and accelerating the chromospheric gas. Locally, in a spicule, the 
energy requirements are comparable to the chromospheric and coronal heating 
requirements; again, it is probably a mistake to separate the spicule problem 
from the overall energy balance of the solar atmosphere. 

In this review we discuss theories which invoke waves to heat the corona. 
If we are willing to interpret the word 'waves’ broadly enough, there are good 
reasons for invoking waves. First, any mechanical process requires that the 
convection zone do work, followed by the mechanical transfer of energy upwards 
into the corona. The solar atmosphere must move rf work xs to be done, and it 
is a fact of life that virtually all motions in the solar atmosphere obey 
hyperbolic equations which yield wave or wave-like solutions. For example, 
the linearized versions of the twisting motions invoked by Parker in the next 
paper obey the Alfvenic wave equation. Second, the corona is obse^ed to 
contain ubiquitous non-thermal motions of the order of 10-30 km s" (rms) (e.g. 
Bonnet, 1978; Cheng et al., 1979; Doschek and Feldman, 1977; Feldman et al., 
1975) . These motions are unresolved in space and time. For the reasons given 
above, it is likely that these motions can be thought of as waves. We will 
argue below that the observed motions may contain sufficient energy to heat 
the corona, and a wave theory of coronal heating seems possible. Third, the 
solar wind may serve as a prototype. Alfven waves (e.g. Belcher and Davis, 
1971) and/or surface waves (Hollweg, 1982a) appear copiously in the solar wind 
beyond 0.3AU and there is some radio evidence for the presence of significant 
wave fluxes in the acceleration region of the wind (Hollweg et al., 1982a). 
Successful Alfven-wave-driven solar wind models have been constructed (e.g. 
Hollweg, 1978a). And the behavior of heavy ions in the solar wind suggests 
the operation of wave-particle interactions (see Isenberg's review in this 
volume) . The solar wind may be telling us that the sun radiates energetically 
significant wave fluxes, and that waves can heat (and accelerate) at least 
part of the solar atmosphere. 

In short, it seems that it should be possible to construct a successful 
wave theory for coronal heating. And in view of the first point in the pre- 
vious paragraph, it is possible that virtually every theory can at some level 
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be thought of as a wave theory. Nonetheless, a successful wave theory has not 
yet emerged. In the following we will summarize some current thinking on the 
subject. We will point out where wave theories succeed and where they fail. 
And we will suggest some possible routes to be followed in the future. 

Some other recent reviews are Hollweg (1981a), Kuperus et al. (1981), 

Priest (1982a), and Wentzel (1978a, 1981). 


Energies Required and Available 


We begin with a brief definition of the problem in terms of the energies 
required to heat the corona. It is useful to split the corona into three 
types of regions: i. coronal holes, out of which high-speed solar wind 
streams flow (e.g. Hundhausen, 1977; Zirker, 1977); ii. quiet corona, con- 
sisting of large-scale closed field regions, such as the helmet streamers ; 
iii. active region loops, consisting of small-scale coronal regions of en- 
hanced pressure with a loop-like morphology, presumably tracing out closed 
magnetic field lines (e.g. Rosner et al. , 1978; Webb, 1981; Withbroe, 1981). 
Convenient summaries of the energy requirements of these regions have been 
given by Withbroe (1976, 1981). 

Coronal holes and quiet corona lose energy via radiation and heat conduc- 
tion back down into the chromosphere. In addition, coronal holes lose energy 
by heat conduction out into the solar wind; this latter energy loss could be 
very large if the high-speed streams are driven thermally, as discussed by 
Olbert in this volume. Ho^es jnd quiet corona require an energy flux density 
of a few times 10 erg cm s entering from below. And if the ^igh-spee^ 
streams are thermally driven, the holes may require as much as 10 erg cm 
s . The volumetric heating rate can be estimated by dividing the energy flux 
density ^by the distance over which the heating occurs. If we take 3x^0 erg 
cm ^ s for^the^flux density, then the heating rate lies between 10 and 

10 erg cm s for heating distances in the range (0.05 - 0.5)r . 

The active region loops lose energy via radiation and via heat conduction 
along the magnetic field lines back down into the chromosphere. In this case 
a useful rule-of-thumb can be obtained by paraphrasing Rosner et al. (1978) 
and Withbroe (1981) . Let be the volumetric heating rate of the 
plasma. ^ The radiation loss out of the (optically thin) plasma is n^ $ erg 
cm s , where n is electron concentration and $ is a function of T 
(electron temperature) only. If we neglect flows, the energy equation is 
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where ^ is heat conduction. Assuming that q is classical electron heat 
conduction along the magnetic field, we can"^ewrite (1) in the form 


1 2 
- 2 % 


K 


E T 
He 


5/2 


dT 



2 h 
p $ T ^ dT 
J *^e e e 


( 2 ) 


where it has been assumed that the loop's cross-gectional area is constant (k 
is Boltzmann's constant, p is pressure, and K T is the heat conductivity). 
We now take $ and p to be constants. The°loop is assumed to have a maxi- 
mum temperature T where q =0, and it is assumed that q -»■ 0 as T -»• 0; 
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The quantity $ can be eliminated in favor of the loop length as follows: From 
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the usual expression for we have 
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where s is distance along the loop. Inserting (4) into (5) and integrating 


from T = 0 to T = T gives 
e e max 
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where L is twice the distance from T = 0 to T 

TT e max 

gives finally 
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Inserting (6) into (3) 


(7) 


Since $ does not appear in (7) , it is probable that errors associated with 
restricting $ to be a constant v/ill not be very significant. Equation (7) 
gives the required energy flux density if all the energy comes up along the 
loop from one of its footpoints; half that value is required if eq^al flu^gs 
come up both footpoints. For a short (long) loop we take L_= 6x10 (6x10 ) 

cm and T =2 x 10° (2.5 x 10°) K, agd wjth = 8.4 x 10 (c.g.s.) we ob- 

tain E L^?^l X 10^ (2.4 X 1^°) erg cm”" s ; th| corresponding volumetric 
heating rates are 1.8 x lo”^ (4 x lo”^) erg cm s . It is interesting to 
note that both the energy flux density and the volumetric heating rate are 
smaller on longer loops; it is the task of theory to explain this. 

Can the observed nonthermal velocities in the corona supply the required 
energy flux densities? It is easy to show that slov; (sound) waves cannot do 
the job, and we will henceforth ignore them, For^the fast or Alfven modes, we 
calculate the energy flux density to be 2p 6v ^ v^ where v^ is the Alfven 
speed and p is the coronal density (the factor^^2rgallows_|or 2 polarization 
states). For coronal holes we take p = 3.3 x 10 gm cm^ , B = B^Cau^s, 

6v = 30 km s” , and the energy flux de^^ity is_^.4 x 10 erg cm. s . For 
active region loops we take p =5.2 xio”"^" gm cm (for a mean^molecular 
weight of h, this density yie?ds a tota^ pressure of 2 dyne cm ^if T = 2^3 x^^ 
10° K) , B = 100 Gauss, 6v = 30 km s , and we obtain 3.7 x 10 erg cm s 
for the energy flux densily? Thes| energy fluxes are adequate to supply the 
required energies, if the 30 km s nonthermal motions are fast or Alfven 
waves . 

Can the nonthermal velocities supply the required volumetric heating 
rates? Write 


E„ = 2m, (2p 6v 
H i o rms 


(8) 


where m. is the imaginary part of the wave (angular) frequency (the second 
factor ^2' again allows for 2 polarization states). Coronal holes require E 
< lO” erg cm”^ s” , and (8) implies m . /w< 0.4 if the wave period is 300s (the 
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latter figure is a guess) . Active region loops require < 2 x 10 erg cm 

s“^, and (8) implies o)./o>< 0.17 if the period is 100s (the reason for this 
choice of period will ^e given below) . Our estimates of u^/co are not larger 
than 1, and it can be meaningful to talk about propagating, but damped, waves. 

Finally, we must ask whether the convection zone can do enough work on the 
system. We will assume that the work is done on the intense photospheric 
magnetic flux tubes (for a review see Spruit, 1981a). If the energy propa- 
gates as an Alfven wave, the convective motions can supply an energy flux 
density to the corona of 


2p 6v ^ V, (B /B ) 
^c c,nns Ac cor c 
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Here the subscript 'c' refers to the top of the convection zone, the factor '2' 
allows for two polarizations, and the factor (B /B ) represents the area ex- 
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pension of the flux tube. Taking p = 3 x 10 ^gm cm rms ^ ^ 

s"^, we obtain energy flux densities of 2*5 x^lO erg cm s lA a coronal 

hole (B =8 Gauss) and 3 x 10 erg cm" s” in an active region loop (B^ 

= 100 Gauls) . These values exceed the requirements by more than an order-ot- 
magnitude. As we shall see, this is fortunate, since most of the wave energy 
is reflected before reaching the corona. (Strictly speaking, equation 9 
ignores the details of the wave generation process at the top of the convec- 
tion zone. A detailed analysis of the coupling between convection and the 
photospheric flux tubes is really needed, but not available.) 

Reflection and Transmission 


Fast waves 


The coronal pressure is small compared to the^magn^ti^ pressure. The fast 
mode dispersion relation is then approximately o) = k v^ , where k is the 

Upon splitting k into horizontal (h) and vertical (v) components. 


wavenumber . 
we obtain 
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Taking the to — k structure for the known solar motions in^the photosphere and 
chromosphere, Ho^lweg (1978b) found that the sun yields k^ < 0 in the corona. 
Fast waves can be expected to be evanescent in the corona, i.e. they suffer 
total internal reflection somewhere below the corona. Fast waves can not be 
expected to supply the required energies to the corona. Moreover, the situa- 
tion is worst in the active regions, which have the greatest energy require- 
ments; the reason is that the active regions are observed to be highly struc- 
tured (implying large k ) and to have strong magnetic fields (implying large 
V ) . Leroy and Schwartz (1982) and Schwartz and Leroy (1982) concur with this 

conclusion. (See also Osterbrock, 1961.) 

However, we note that interesting results have been obtained by Habbal et 
al. (1979). They postulate the presence of coronal fast waves with a period 
of 3s. They use ray-tracing techniques to follow the propagation of the wave 
energy, and emphasize the tendency of the fast waves to refract into regions 
with smaller v ; some waves cem even be trapped inside a dense loop, in 
analogy with tfte trapping of light in an optical fiber. They consider the 
Landau/transit-time damping of the waves (e.g. Barnes, 1966). ^ An interesting 

feature of the damping is that it increases with 6 = 8irp /B , if 3 is small 
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(the subscript *p* refers to the protons). In the isothermal model of Habbal 
et al., 3 is largest where v is smallest. Thus the waves refract toward 
regions o? larger damping. This means that the fast waves can heat some 
coronal regions more strongly than others. The effect is enhanced by a 
'positive feedback* , whereby local heating increases 3 along a given field 
line (since heat conduction spreads the heat along B) ,^which in turn increases 
the heating, and so on. A 'catastrophic' situation can occur in which once 
fast waves begin heating the plasma, they dump all their energy in a small 
flux tube. Habbal et al. suggest that such a scenario can account for the 
highly structured nature of the corona in a natural way. But it has to be 
shown that the waves exist, and that equation (10) can be overcome. 

Alfven waves 


Unlike fast waves, these waves never totally internally reflect, and we 
consider them in some detail. Hollweg (1981b) has considered the propagation 
of small-amplitude (linearized) axisymmetric twists on a background potential 
magnetic field which has an axisymmetric untwisted fleur-de-lis structure. If 
the axis of symmetry is vertical, the twisting motions do not couple to 
gravity or to the radiation field, and they are non-compressive . If the 
motions vary as exp(io)t), they obey 
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where x = 6v^/r, y = r the subscript 'o' refers to the background, the 

prefix '6* refers to the wave, 0 is the angle about the axis of symmetry, s is 

distance along any field line, and r(s) is the distanc^ from the axis to the 

field line. For field lines near the axis we expect r B 'X/ constant, and (11) 

, o 'v 
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Equation (14) is the Alfvenic wave equation; we deal with it in what follows. 

If V 'X/ e , then (14) has the solution 
/\ 
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where ^ = 2hio/v^(s), and a and b are complex constants. From (13) 
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<S>r is along the magnetic field: 
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From (15) and (16) we obtain 
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From the form of ,118), we interpret the part of (15) as the upward-going 

wave, and the H ^ part of (15) as the downward-going wave (Hollweg, 1972). 

Now the exponential behavior of v used in (15) - (18) roughly represents 
the probable behavior in the chromosphere and transition region of a flux 
tube, but V is much more nearly constant in the coronal part of a flux tube 
(Figs. 3 an§ 10 of Hollweg, 1981b). If the corona extends to infinity, as in 
a coronal hole, it is possible to obtain some useful analytical results by 
considering a two-layer model, in which 


V = V = constant, s > 0 (19a) 

A A , cor 


V, = v^ e 
A A, cor 
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Here s > 0 represents the corona, and s<0 represents the chromosphere and 
transition region. Equations (15) - (18) apply in s < 0, while the usual 
harmonic solutions apply in s > 0. At s = 0 we apply the matching conditions 
that 6v and 68^ be continuous (this ensures continuity of <S^>) . The 
matching conditions give 
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where the argument of the Hankel functions is 2ho)/v = a. From (20) we 

can calculate the reflection coefficient, R, i.e. tn^ ratio of downward-going 
energy flux to upward-going: 

R=lbl2/la|2 (21) 


In the limit of small a, we obtain 

R l-2iTa 

(Equations 19-22 have been given by Leer et al. (1982) . They are a special 
case of Hollweg (1972).) Since equations (20) - (22) are obtained from the 
matching conditions at s = 0 it seems natural to think of the reflection as 


10 


occurring at the discontinuity in scale height there (Hollweg, 1972) . 

However, Wentzel (1978b) , following Tolstoy (1973) , and Leroy (1980) have 
argued that R is really the accxmiulated reflection off of the entire chromo- 
sphere and transition region. The latter interpretation is probably correct, 
but we are then perplexed as to why R can be calculated solely from the 
matching conditions at s = 0. 

For a flux tube_|ntering a coronal hole we might have h = 250 km, and 

cor ~ 1200 km s . If the wave period is 300s, the energy transmission 
coefficient is 1-R = 0.055. In the previous section we estimated that the 
convective work c|n_|upply 2.5 x 10 ' erg cm s” . Multiplying by 1-R gives 
1.4 X 10 erg cm s available to the corona. This is enough to supply the 
energy requirements. 

The situation is different on an active region loop. We imagine the loop 
to be stretched out so that (possibly important) complications arising from 
the curvature are ignored (e.g. Wentzel, 1978b) . But we must now at the very 
least consider a 3-layer model, since there is a chromosphere and transition 
region at each end of the loop (an 18-layer model was considered by Hollweg, 
1981b) . We take 


= '^A,cor ^^P<s/2h), s < 0 


(23a) 


V, - f 0 < s < L 

A A, cor 


^A = '^A,cor r-(s-L)/2h'], s > L 


(23b) 

(23c) 


We imagine there to be a source at the far end of the loop, in s > L. We 
impose the boundary condition that there be only an outgoing (from the source) 
wave in the region s < 0. The other two regions have both upward-and 
downward-going solutions. Applying matching conditions at s = 0 and s = L 
allows us to calculate the ratio of downward-going energy flux to upward-going 
flux in the region s > L. Again calling this ratio R, we obtain 


R =|N/Dj 
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We have defined ot= 2ho)/v, , 0= 2h* w/v, , and k = 0 )/v, 

A, cor A, cor A, cor 
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Expressions (24) - (26) have not yet been analyzed in detail. But they 
take a simpler form if a = 6, and if o, 3 << 1. Then 

N % Y,^sin kL + 2 Y Y, cos kL (27) 

^1 o 1 


and 


D N + 2iY^ cos kL (28) 

It is interesting to note that N can be exactly zero at certain resonant 
frequencies which are approxmately given by 
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where m = 1,.2, 3... (For example, the resonant period is 40s if m = 1, 

= 3000 km s” , and L = 6 x 10^ )on.) At these frequencies the reflections ' 
vanish, and a large energy flux can pass through the corona and out the other 
end of the loop. In the vicinity of one of the resonances it is possible to 
use (24) , (27) and (28) to obtain the following expression for the energy 
transmission coefficient: 
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Denoting the full-width-at-half -maximum by Aw, we find 
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(Q denotes the quality) . If L = 6 x 10 km and h = 150 km, then Q = 32; the 
resonance is moderately high quality. 

We have already estimated that the convection zone work can supply 3 x 10 
erg cm s . If this power is in a band width B , then the energy flux 
density passing through the corona in one resonant peak is 


F 

res 


(3 X 10®) 


Aw 



ir 

2 


(32) 


where the factor (tr/2) comes from integrating the area under (30). From Fig. 
3 of lonson (1982) , we estimate B 3 x 10 ^s . Then from (29) and 
(31) -(32) we find “ 
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Note that longer loops receive less pjwer, as (qualitatively) observed. If we 
take h = 150 “2^-2^ km s , L = 6 x 10 km and m = 1, then 

= 7.8 X 10 4rg cm s . This is more than enough energy to supply the 
loop, but this energy goes right through the loop, since we have put no damp- 
ing into the calculation. The calculation with damping still needs to be 
done, but work in progress shows that adequate energies can reach and stay in 
the loop, since damping tends to broaden the resonance peak, compensating for 
its reduction in height. 

Equations (25) -(33) are new, but some of these ideas have been considered 
previously. lonson (1978) was the first to show that standing waves can be 
excited on active region loops. The first paper to show that resonances can 
eliminate the reflections and allow Alfven wave energy to enter the corona was 
by Hollweg (1981b). The same point v;as made independently by Zugzda and 
Locans (1982) , but their Alfven wave equation differs from (11) and (12) 
above. lonson (1982) has discussed the resonances using RLC circuits as an 
analogy. This approach omits some physics, however. For example, ^in the 
absence of dissipation lonson gets Q = ", and he uses Q = 10'^ - 10^ in his 
paper; by contrast, we find Q = 30 or so, even in the absence of dissipation. 

An interesting point made by lonson (1982) is that there appear to be two 
classes of loops (active region loops and * large scale structures'). There 
may also be two peaks in the photospheric power spectrum. lonson suggests 
that a narrow-band resonance is a natural means by which the details of the 
power spectrum could be mapped into the corona, and thus a double-peaked 
spectrum can give rise to two classes of loops. 

Alfvenic Surface Waves 


The corona is highly structured. In the limit that the structuring takes 
the form of discontinuities, it is possible to find new wave modes supported 
by the surfaces. These surface waves have recently been of considerable 
interest. See Edwin and Roberts (1982) , Gordon and Hollweg (1983) , Hollweg 
(1982a), lonson (1978), Roberts (1981a, b) , Wentzel (1979). 

Consider a cold, stationary background plasma which varies only in the x “ 
direction. The background magnetic field vector varies in its magni- 

tude is constant. The system supports small-amplitude fluctuations obeying 
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the fluctuations vary as 


- and e = 4 ttp^(o^ - 

y + ik z - io)t] . Now 2 assume that everything is uniform except at a 
discontinuity at x = 0* if ^ ^0 everywhere, it is possible to find solu- 
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tions which evanesce away from the discontinuity as expt±qx], 
relation for these surface waves follows by requiring that 6v 
(e/q^)96v /3x be continuous at x = 0« * 

For exctmple, suppose B ^s constant and in the z-direction. Suppose also 
,. 1 . ■_ , __ ^ dispersion relation turns out to be 


that k is large, so that 
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where p is the average of the densities on the two sides of the discontinu- 
ity. Equation (35) is similar to the dispersion relation for Alfven waves. 
These surface waves generally turn out to be rather similar to Alfven waves. 
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In particular, they propagate energy along B . Much of what was said about 
Alfven waves in the previous section applies to surface waves as well. They 
may be suitable candidates for coronal heating, although many details, such as 
their ability to propagate energy through the chromosphere and transition 
region, remain to be worked out. 

It is interesting to note that surface waves may be present in the solar 
wind (Hollweg, 1982a). 

Wiggles of Thin Flux Tubes . 

Spruit (1981a, b) has considered transversal oscillations of thin vertical 
magnetic flux tubes imbedded in a field-free gas. If the tube is in pressure 
and temperature equilibrium with its surroundings, then the (small) horizontal 
displacements of the tube obey: 
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where g is the gravitational acceleration, and 6^ = constant is the ratio of 
gas to magnetic pressures inside the flux tube. If the atmosphere is 
isothermal, we can take ^'v^exp [ikz-iwt] and (36) yields: 


0) ^ 
c 

where 

2 2 

= --g P 

8p (23^+1) 


(37) 


(38) 


4 

The waves have a low-frequency cutoff. If the temperature is 10 K, if the 
molecular weight is 1.3, and if = If we find from (38) that periods longer 
than 900s are evanescent; this is not a severe restriction. 

These waves carry energy along the magnetic field. If oo >> o) , the phase 
and group velocities are comparable to the sound speed. They could carry 
energy into the low chromosphere. However, the tubes cease to be thin above 
the low chromosphere, and a more detailed analysis is necessary. They could 
couple some energy into the fast and Alfven modes, which could carry the 
energy to greater heights. 


Dissipation 

The general conclusion of the previous section is that Alfven or Alfvenic 
surface waves, and to some extent transversal tube waves, can in principle 
supply the corona with its energy requirements. This is the first 
requirement of a coronal heating theory. The second requirement is that the 
waves deposit energy as heat. Wave theories of coronal heating have so far 
failed in this regard, but we will summarize some of the possibilities. 

One point to note at the outset is that the observed coronal motions, if 

interpreted as magnetic waves, are quite linear. In a cgronal ho^e we might 

have v^ 1200 km s (if B = 8 Gauss and p = 3.3 x 10^ gm cm ) , while in 

an active region loop we might have v 4000 km s" (if B = 100 Gauss and p = 

A 'v 
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In all cases 6v /v^ << 1, in contrast to the solar 
riPs A 


Viscous Heating. 


Coronal viscosity is mainly due to the protons. Following Braginskii 
(1965) , the proton collision time in an electron-proton plasma is 


T = 0.75T n 
P P P 


(39) 


for a Coulomb logarithm of 22. In the corona o) x >>1 (w is the proton 
cyclotron frequency) and the protons are well-tiid^to the^field lines. Of the 
five viscosity coefficients given by Braginslcii, n 10 (cgs) is by 

far the largest; shear viscosity is smaller by (w )” (roughly). If only 
Ti^ is considered, the viscous heating rate is ^ 
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where^z is along the magnetic field. For purely parallel motion 

(V*v) while for purely transverse motions Q = r\ (V • v)^ /3. 

thumb we will take Q=n(V*v). ^ ^ 

.P o 'v -3 -3 -1 

In an active regi6n loop we observe Q < 2 x 10 erg cm s 

equation of continuity, this requires 6p^ /p <0,8 if T =2.3 

the wave period is 100s. In a coronal hoTI we observe Q 10~ 

This requires <5p^ /p < 0.9 if T = 1.5 x 10^ K and if ?he wave 

300s. Such large^lensity fluctuations may occur occasionally in 

in the corona, but they are not compatible with the requirements 

the entire corona by any of the waves discussed so far. 
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Heat Conduction 


Coronal heat conduction is due mainly to the electrons. They too are tied 
to the field lines and the heat conduction is along the magnetic field. The 
plasma heating rate due to heat conduction damping of waves is 
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-7 5/2 

where the heat conductivity isic =8,4x10 T (cgs) for a Coulomb 

logarithm of 22. ® 

Low-frequency waves induce nearly adiabatic temperature fluctuations, i.e. 

1 5T I ^ (y “ 1) T |v • 6;^l/o), where y = 5/3 is the ratio of specific heats. 

We then find (takfna T = T ) 

' e p 

^ = 8.4 X 10^ (y-D^T^k^^/o)^ ^^2) 

0 

If = 2 X 10 Kj we find that viscous heating dominates heat conduction if 
t»/k^ > 860 km s . Coronal Alfven speeds generally exceed this value, and we 
deduce that heat conduction will be even less effective than viscosity in 
damping the waves. 

For completeness we should mention that heat conduction smooths out the 
temperature fluctuations in high-frequency waves. This reduces the 
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effectiveness of heat-conduction dampinq even further. See Gordon and Hollweg 
(1983) . 


Shocks. 


We have already mentioned that Alfven, fast, or Alfvenic surface waves are 
to a good approximation linear in the corona. Shocks probably do not form in 
the corona. However, Hollweg et al. (1982b) have suggested the possibility 
that shocks can form in the chromosphere, on their way to the corona. Hollweg 
et al. considered Alfven waves, which steepen into a train of shocks which are 
nearly switch-on shocks (see, for example, Boyd and Sanderson, 1969). These 
shocks enter the corona, and can carry substantial energy fluxes. Hollweg et 
al. suggested that some of the impulsive events observed in the transition 
region (see the article by Dere in this volume) could in fact be the shocks. 

The volumetric heating rate due to a periodic train of weak switch-on 
shocks with period t is 

B 2 

Q = 

SOS 32 ttt 

if the coronal pressure is small compared to the magnetic pressure ( Hollweg, 
1982b) ; Av^ is the j^p in the velocity component transverse to the shock 
normal. Since Q , this mechanism yields negligibly small heating rates 

in active regions. However, it could conceivably work in coronal holes or 
quiet corona if and t are small enough, and if Av is large enough. For 

example, if = 5 Gauss^^p = 3.3^^ 10^ gm cm , Av = 200 km s , and t = 
100s, we obtain ^0 €^rg cm s . The corona can be heated over an 

extended distance by this mechanism. 

Switch— on shocks are probably a worst case. A best case for shocks is a 
train of fast shocks propagating across the magnetic field. If the shocks are 
weak and if the coronal pressure is small compared to the magnetic pressure, 
the volumetric heating rate is 


Av 4 

( ^ ) (43) 
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where Av is the velocity jump across the shocpk. In an active region loog we 
might have B^=100 Gauss and v^^ =.|000 km s^*". If we take Av = 200 km s ^ 
and T = 100s we find Q = 6 x^lO erg cm s . This is comparable to the 
heating requirements of some moderate length loops. In a coronal hole we 

^ ^ Gauss and v^^ = 800 km s“^, and we find Q = 7.8 x lo"^ erg 
cm s . The heating of the coronal hole is substantial in fhis case. 

Thus fast shocks propagating across B can yield substantial heating. The 
same is presumably true for fast shocks^prcpagating at some not-too-small 
angle to B^. But where would such shocks come from? One guess is that they 
could form in the chromosphere in the manner investigated by Hollweg et al. 
(1982b) for switch-on- shocks, in fact, the study of Hollweg et al. is 
probably unrealistically restricted, since there is no reason to expect that 
the sun will yield shock normals which are nearly aligned along the magnetic 
field; in view of the strong cross-field structuring observed on the sun, the 
opposite is probably the case, and further studies must be done. In any 
event, we will ultimately have to rely on observations to tell us whether 
shocks with the required frequency and amplitude do in fact form in the 
chromosphere and enter the corona. 
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Phase Mixing and Turbulence . 


Consider two neighboring magnetic field lines along which energy is propa- 
gating in the form of an Alfven wave. If the phase velocities are different 
on the two field lines, then the motions on those field lines will not always 
maintain the same phase relationship. The motions will move in and out of 
phase, and (in the latter case) extremely large cross-field velocity gradients 
can develop. Heyvaerts and Priest (1982) (see also Priest, 1982b) have 
pointed out that this phase mixing is ripe for viscous dissipation. They 
consider the special case where B is in the z- direction, 6v is in the 
y-direction, and v^ = v^(x)- UnfSrtunately, since the motions in this case 
are shears, the viscous damping will be weak, because o) x >>1; this limita- 
tion of the efficacy of viscosity was not considered by'^fiey^aerts and Priest. 

Heyvaerts and Priest consider also the possibility that large velocity 
shears between neighboring field lines can result in Kelvin-Helmholz instabil- 
ities. (They consider in particular the situation on coronal active region 
loops, where the resonances discussed in section IIIB lead to nearly-standing 
waves on the loops. A similar idea was mentioned earlier by Hollweg (1981b) 
with regard to Alfven waves in the chromosphere. Since the solar atmosphere 
IS not homogeneous, and since there is no reason to expect phase coherence 
between motions on different field lines, phase mixing may be a ubiquitous and 
important effect.) Priest (1982b) estimates the Kelvin-Helmholz growth rate 
to be 0 )^ kj^|6y|, v;here is the transverse wave number associated with the 

velocity shears. But it must be kept in mind that the driving velocity shears 
in this problem are time-dependent. For the analysis (which assumes that the 
driving shears are steady) to hold, it is necessary that w. >> w. If we take 
0 ) = 2tt/(100s), and |6y| = 40 km s (corresponding to 6v ^ 30 km s~^) , we 

require = 2ir/kj^<< 4000 km. This is not a strong consirlint (the active 
region loops are only a few thousand km in diameter) , and Kelvin-Helmholz 
instabilities may occur. The instabilities may initiate a turbulent cascade 
to higher wavenumbers where viscosity (or some other process) can convert the 
energy into heat. 

Heyvaerts and Priest (1982) and Priest (1982b) consider also the possibil- 
ity that the phase relationship between Alfven waves on neighboring field 
lines is such as to produce large magnetic shears. Tearing-mode instabilities 
are then possible. (This idea was also discussed earlier by Hollweg, 1981b.) 
Priest (1982b) estimates the tearing growth time to be 


tmi 


% |«yI 
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n k, 


(45) 


where n 2 x 10 T^ (cgs) is the m.agnetic diffusivity; in computing g we 
have taken the perpendicular electrical conductivity (Braginskii, 1965) and a 
Coulomb logarithm of 22. (Note that ri is so small in the corona that classi- 
cal electrical resistivity fails miserably as a dissipation mechanism. ) For 
the analysis (which assumes that the magnetic shears are steady) to be valid, 
it is necessary that wt « 1. If the wave period is 100s, this requires X 
« 10 km, where we have taken T = 2.3 x 10 K, and |6y| = 40 km s"^. It is ^ 
not known whether the coronal waves are structured on these scales. But 
tearing instabilities, if they occur, can initiate a turbulent cascade to 
higher wavenximbers. 

Suppose a turbulent cascade is initiated. How effective will it be? 
Unfortunately , most turbulence theory has been developed for incompressible 
fluids with isotropic turbulence. Neither condition can be expected to apply 
in the corona. However, as Montgomery (this volume) says, "it's the only game 
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in iown” , so we'll play anyway. For fully-developed Kolmogoroff turbulence 
one expects 

E(k) (46) 

where k is wavenumber, % 1^5 is a_^niversal (?) constant, /E(k)dk = <6v^>, 
and e (dimensions = velocity time ) is the rate at which <6v^> cascades to 
high wavenumbers. Integrati,ng (46) from k^ to « gives 
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2<6v^> 


k <5v^> 
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2(1.5 C^j)3/2 


(47) 


where it has been as.sumed that all the energy is in k > k . Now we can regard 
e/( 2 < 6 ^ 2 >) ag ^ measure of the damping rate, w. , since the cascade pumps 
energy to arbitrarily high wavenumbers where it is absorbed and converted into 
heat. Equation (47) then yields (c.f. Section II) 
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(we have again allowed for 2 polarization states) . The problem is that we 
are not sure what to put in for k . If we take k = with 3000 km (a 

reasonable guess^for an active region loop), p = 5.2 x 10 gm cm and 

^^rms “ ^ , we find = 2.5 x 10 erg cm s . Comparison with 

Section II shows that this is adequate to heat the active region loops. In a 
coronal hole we mighty ha^e to be the mean distance between the photospheric 
magnetic fl^x tubes^^'ilO^ km, and with p = 3.3 x 10 gm cm we obtain E 
= 4.7 X 10 erg cm s . Again, we obtain an adequate heating rate. We ^ 
tentatively conclude that a turbulent cascade can provide the required heat- 
ing. But a theory for turbulence which is applicable to the corona must be 
developed. We suggest that this subject be vigorously pursued in the future. 
(The author thanks Dr. C. Smith for advice on this paragraph.) 


Surface Waves 


lonson (1978) has considered the propagation of surface waves on 
non-discontinuous 'surfaces*. He considers the problem discussed previously, 
which has (35) as its dispersion relation in the case of a truly discontinuous 
surface. He finds, however, that m is now complex, with 



where *a' is the thickness of the 'surface* and Av is the difference between 
the Alfven speeds on the two sides. lonson interpreted o). as a damping rate, 
but Lee (1980) pointed out that it is not correct to thin^ of damping of a 
normal mode, because there are no normal modes in this problem (unless a = 0) . 
Instead, the appearance of represents a readjustment of the system's energy 
distribution. (Lee suggests that the situation is analogous to leakage of 
particles out of a potential well, as in radioactive decay, (o. represents the 
leakage out of the well, but there is no net loss of particles^) The energy 
flows into the surface. Steep gradients develop in the surface, and viscous 
dissipation will eventually occur in a thin layer. It seems likely that the 
dissipation in the thin layer will adjust itself to absorb the energy flow 
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into the surface, as suggested by lonson (1978). The problem with dissipation 
still needs detailed study, but it may represent a promising mechanism for 
converting wave energy into heat in many thin layers in the corona. (See also 
Lee et al., 1983). 

Summary 


The corona is observed to move with velocities of 30 km s ^ or so, rms. 
Since most coronal motions obey hyperbolic equations, it is reasonable to 
think of the motions as being waves. 

If the motions are Alfven or Alfvenic surface waves, they can carry the 
required energies into the corona. Some large fraction of the energy is 
presumably reflected in the chromosphere and transition region, but theory 
indicates that sufficient energy can be transmitted into the corona. Loop 
resonances may play a special role in allowing energy to enter the active 
region corona. 

The problem is with the dissipation. Viscosity is the most promising 
dissipation mechanism, but special conditions have to be fulfilled. The wave 
energy must ultimately appear at sufficiently short spatial scales so that 
viscosity becomes effective. Fast shocks may be suitable, since the shock 
thickness automatically adjusts itself to yield the required entropy jump 
across the shock. A turbulent cascade seems even more promising, but an 
anisotropic compressive turbulence theory still needs to be formulated. 

Surface waves may be yet another route by which the wave energy finds itself 
at small enough spatial scales for viscous heating to be effective. But the 
detailed analyses of these processes still need to be done. 

Another area where more work is needed is the following: We have conceptu- 

ally separated the coronal heating from the chromospheric heating. And we 
have regarded the chromosphere and transition region as being fixed entities 
which carry (and reflect) the waves. But the coupling between the chromo- 
sphere, corona, and transition region may be an integral part of the entire 
atmosphere's energy balance. And the coupling may be strong and dynamic. For 
example, Hollweg (1981b) and Hollweg et al. (1982b) have argued that the 
chromosphere and transition region may themselves be set into vigorous motion 
as the waves propagate from the photosphere into the corona. (A possible 
connection with the spicules has been noted by Hollweg et al. (1982b) and by 
Hollweg (1982c).) It is not yet clear how this dynamic coupling affects the 
ideas presented in this review. 

Finally, we have been implicitly regarding the waves as non- impulsive. Yet 
the observations discussed by Dere in this volume suggest the opposite: the 

transition region moves violently and impulsively. What are the implications 
of these observations for coronal heating? We have already mentioned that the 
impulses could be shocks, but the issue is far from being resolved. These 
observations demand further study. 
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